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Abstract 

We have examined the effect of a c(2x2) overlayer of CO on the surface phonons of the substrate, 
Cu(OOl), by applying the density functional perturbation theory with both the local (LDA) and 
the generalized-gradient (GGA) density approximations, through the Hedin-Lundqvist and the 
Perdew-Burke-Ernzerhof functionals, respectively. Our results (GGA) trace the Rayleigh wave 
softening detected by helium atom scattering (HAS) experiments to changes in the force constants 
between the substrate surface atoms brought about by CO chemisorption, resolving an ongoing 
debate on the subject. The calculated surface phonon dispersion curves document the changes in 
the polarization of some modes and show those of the modes originally along the Y direction of the 
clean surface Brillouin zone (SBZ) which are back-folded along the A direction of the chemisorbed 
SBZ, to be particularly consequential. The vertical and shear horizontal section of in the SBZ 
of the clean surface, for example, is back-folded as a longitudinal-vertical mode, indicating thereby 
that Si — predicted a long time back along A for the clean surface — may be indirectly assessed at 
X upon CO adsorption by standard planar scattering techniques. These findings further suggest 
that some of the energy losses detected by HAS along A, which were associated to multiphonon 
excitations of the adlayer frustrated translation mode, may actually correspond to the back-folded 
substrate surface modes. 

PACS numbers: 68.35.Ja, 63.20.D-, 68.43.-h, 63.22.-m 
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I. INTRODUCTION 



Although it is understandable that a large part of theoretical investigations in catalytic 
surface science is dedicated to developing an understanding of the surface electronic and 
geometric structure and energetics of processes such as chemisorption and thermal activa- 
tion, as a function of catalyst element and surface geometry, two decisive aspects are often 
overlooked: attestation of the dynamical stability of the model system^ and understand- 
ing of the vibrational dynamics of the react ant- catalyst complex. Clearly, analyses lacking 
considerations of system dynamics disregard the fact that reaction paths and the so-called 
attempt frequencies or rate preexponential factors pertaining to dynamical processes are 
themselves determined by the displacement patterns and frequencies of the phonons of the 
system under consideration, respectively. Moreover, according to the harmonic transition- 
state theory, the entire spectrum of phonons (at the equilibrium and the transition states) of 
the composite system is required to determine the attempt frequency of any given process, 
and not just those vibrational modes in which reactants are primarily involved. Regard- 
less, it is frequently assumed that the changes in the free energy, which govern the rates of 
adsorbate processes, involve little contribution from the substrate because of relatively low 
chemisorption energies, as in the case of carbon- monoxide (CO) adsorbed on noble metal 
surfaces. The CO binding energy on Cu(OOl), for example, is 20 times smaller than the 
carbon-oxygen binding energy and about three times smaller than the chemisorption en- 
ergy of CO on Ni(lll).-i^ Nevertheless, examination of ultra-violet and X-ray photoelectron 
spectra signalize a significant interaction between the molecular and the Cu orbit als.-'^'^'^'^ 
Namely, the valence levels of CO molecules adsorbed on Cu are rearranged, broadened, and 
shifted with respect to those of molecules in the gas phase, in a similar fashion to that 
observed for Ni(OOl) — albeit, of course, to a lesser degree. In each of these cases, molec- 
ular adsorption also impacts the substrate phonons. One such example is the case of the 
hydrogen overlayer on Pt(lll),^°'^^ in which the frequency of the substrate Rayleigh wave at 
the zone boundaries is modified substantially from its value on clean Pt(lll). Furthermore, 
even in the case of Cu self-diffusion,— in which the substrate may be expected to play a 
less important role due to the short-range interaction among coinage metal atoms, it has 
been shown that the contribution of the substrate to free energy changes in the course of 
adatom hopping is non negligible. Should the Cu-CO interaction have a range comparable 
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to or longer than the Cu-Cu interaction, for instance, grasping the microscopic details of 
CO-related reactions facilitated by metallic surfaces will demand taking into consideration 
all inter-atomic couplings within the system. 

We have recently applied the density-functional perturbation theory^^ii^ (DFPT) to ex- 
amine the vibrational modes of the c(2x2)-CO adlayer on the Cu(OOl) surface. Since all 
pertinent interatomic interactions are automatically included in our approach, we were able 
to calculate the phonon frequencies at arbitrary propagation directions, i.e., in the en- 
tire surface Brillouin zone (SBZ).— Analysis of the displacement vectors showed that the 
CO-modes are influenced by molecule-substrate and molecule-molecule interactions.— Inter- 
actions among neighboring CO molecules, separated by ~ 3.6A, were such as to disperse 
and/or split the C-0 stretch, the frustrated rotation (FR), and the frustrated translation 
(FT) modes. Interestingly, we also found that omission of the dynamics of the substrate in 
the calculations lowers the frequency of the Cu-CO stretch mode by ~8 meV with respect to 
the value obtained from DFPT calculations. The frequencies of two adlayer modes, the Cu- 
CO stretch and the FT modes, were found to depend strongly on whether the local density 
(LDA) or the generalized gradient (GGA) approximation, as formulated by Perdew, Burke, 
and Ernzerhof (PBE), was used. In fact, the results are a testimony to the unsuitability of 
LDA to describe the CO adlayer since, contrary to experimental assessments,->i^ not only 
does it render the top adsorption site as a shallow local minimum,- but also implies that 
the FT mode of the adsorbed CO molecules is unstable almost everywhere along the A 
and S directions of the c(2x2) SBZ (Fig. [1] (b)). Such dynamical instabilities are indicators 
of the inability of LDA to describe the Cu-CO interaction which consequently leads to er- 
rors in the predicted CO adsorption site on the surface. The problem, in turn, originates 
from the expression for the exchange-correlation energy, which is inherently approximated 
in the Kohn-Sham formulation and thus leads to a non-exact cancelation of the Coulomb 
self-interaction.— GGA-PBE, though not systematically self-interaction free, introduces an 
enlargement of the energy gap between the highest occupied orbital (HOMO) — 5a — and 
the lowest unoccupied orbital (LUMO) — 27r* — of CO, thus reducing the hybridization 
between the 27r* orbital and the metallic (i-states,— removing the discrepancy between 
theory and experiment regarding the preferred adsorption site of CO on Cu(OOl),^ and ac- 
curately reproducing the dispersion of the FT mode as measured by helium atom scattering 
(HAS) experiments.— 
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In the present work, we turn again to the real-space force constants and the phonon 
dispersion of c(2x2)-CO/Cu(001) to analyze the effect of CO molecules on the dynamics of 
the Cu(OOl) substrate. First of all, the different scenario exhibited by GGA-PBE regarding 
the C-Cu interaction calls for revisiting the effect of the CO adlayer on Cu(OOl). Ellis et al.— 
observed via HAS measurements that the Rayleigh wave (RW) of Cu(OOl) softens upon CO 
adsorption with respect to the clean surface. The later effect was only partially explained 
by the mass overloading of the CO- covered Cu atoms since the softening obtained by simply 
increasing the mass of such atoms while keeping intact the force constants of the clean surface 
underestimates the observed softening. Since no significant changes in the force constants 
of the Cu substrate was found in the " frozen-phonon" LDA DFT calculations of Lewis and 
Rappe^^ who obtained reasonable agreement with the experiment for the frequency of the 
back-folded RW at F, the mass overloading effect was accepted as the main reason for the 
RW softening. In reality, it is not straightforward to interpret or classify the softening as 
due to either factor. On the one hand, while only one of the two Cu surface atoms in the 
unit cell adsorbs and carries CO, the RW refers (mainly) to the vibration of the first layer 
which is represented by both atoms in the surface unit cell. On the other hand, covered and 
bare atoms relax in opposite directions upon CO chemisorption. In other words, covered 
atoms not only support a CO molecule but also relax outwards, whereas bare atoms relax 
further inwards with respect to bulk interlayer spacing (GGA-PBE). Developing a rationale 
for the softening of the RW becomes even more complex given that despite the disagreement 
between LDA and GGA-PBE in the predicted structural features and changes in the force 
constants of the substrate — , both reproduce reasonably well the HAS measurements of the 
RW at the zone center. 

In this work we also focus on the substrate modes which exist along the Y direction of 
the clean SBZ and are now back-folded along the A direction of the chemisorbed SBZ with 
changed polarization. Our calculations in fact suggest that some of the energy losses detected 
by HAS in this direction — and ascribed to the multiphonon excitations of the adlayer FT 
mode REF — may actually correspond to surface back-folded modes. Such backfoldings 
and polarization changes may be the key for the experimental detection of modes that 
are otherwise inaccessible to planar scattering techniques. Of particular relevance is the 
vertical (V - vibration perpendicular to the surface) and shear horizontal (SH - vibration 
perpendicular to the propagation direction) section of Si predicted for the clean surface. 
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which may unfold 5*1 at X to planar scattering detection. Specifically, such a mode is back- 
folded as a longitudinal (L - vibration pa rallel to the propagation direction) and vertical 
mode for the chemisorbed surface and is degenerate with the pure shear-horizontal mode, 
Si, at the zone boundary. 

The rest of this work is organized as follows: Section 11 contains the computational details. 
Section 111 is a summary of results concerning the structure of bulk Cu, the clean Cu(OOl) 
surface, and the c(2x2)-CO/Cu(100) chemisorbed surface. In Section IV, we present our 
results and discussion of the dynamics of all three systems. Finally, Section V contains 
concluding remarks of this study. 



II. COMPUTATIONAL DETAILS 



Periodic super-cell calculations are performed on the basis of the DFT formalism and 
the norm-conserving pseudopotential approach.— The present results are derived from the 
mixed basis (MB) technique.— Results using both LDA and GGA are obtained. The for- 
mer is applied through the Hedin-Lundqvist^ parameterization of the exchange-correlation 
functional, whereas GGA is introduced via the PBE expression.— 

The clean and the CO-chemisorbed Cu(OOl) surfaces are simulated with symmetric slabs 
inside a tetragonal unit cell containing either 9 (for LDA) or 7 (for GGA-PBE) layers of 
Cu. On the chemisorbed surface, CO molecules are symmetrically located on each side of 
the slab. Periodically repeated slabs are separated by a distance equivalent to 11 and 9 
vacuum layers, correspondingly. Integrations inside the Brillouin zone are performed over a 
discrete mesh of 8x8x8 k-points for bulk Cu and of 8x8x1, and 6x6x1 k-points for Cu(OOl) 
and CO-c(2x2)-Cu(001), respectively. 

The calculation of the lattice dynamical matrices at specific q-points of the SBZ is based 
inear response theory embodied within DFPT,— as implemented in the MB scheme 



on the 
in Ref. 



24J . The dynamical matrices for bulk Cu, Cu(OOl), and c(2x2)-CO/Cu(001) are 



calculated at the q-points of a 4x4x4, a 4x4x1, and a 2x2x1 mesh, respectively. The real- 
space force constants in these systems are obtained by the standard Fourier transform of the 
corresponding dynamical matrices.— The force constants of both surfaces are then matched 
with those of bulk Cu to model a clean and a chemisorbed asymmetric slab of 50 layers and 
used to obtain the frequencies at arbitrary q-points. Surface modes on the clean surface 
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have been identified as those whose amphtude weight in the two outermost layers is larger 
than 20%. On the chemisorbed slab, surface modes and resonances have been identified as 
those whose amplitude weight in the 6 outermost atoms (including C and O) is larger than 
20 and 5%, respectively. For further details of the computational methods, we refer the 



reader to Ref. 



20 



III. RESULTS AND DISCUSSION OF STRUCTURAL PROPERTIES 

Our results for the bulk Cu lattice parameter (3.57 A [IDA] and 3.68 A [GGA-PBE]) 
are in good agreement with all-electron^ (AE) and pseudopotential^»^ calculations. Never- 
theless, our calculated bulk modulus, B (170 GPa [IDA] and 128 GPa [GGA-PBE]), falls 
below that provided by AE calculations.^^ Discrepancies in this respect are, in fact, expected 
since B involves the second derivative of the energy with respect to the volume, being thus 
more susceptible to the differences between AE and PP calculations than the lattice param- 
eter. As for agreement with experiment,— i^^ii^ LDA underestimates the lattice constant and 
yields larger bulk modulus. GGA-PBE overcorrects LDA, though it reproduces better the 
experimental bulk modulus than does LDA. Further details and comparisons can be found 
in Ref, 20. 



A schematic top view of the structure of Cu(OOl) and the c(2x2)-CO adlayer on Cu(OOl) 
is shown in Fig. [I](a). The relaxation of the interlayer distances normal to the surface of 
the clean Cu(OOl) has been extensively studied theoretically.-'^'^^'^'^ Nevertheless, before 
conducting our study on the CO- chemisorbed surface, it is essential to test the applied 
methodology on the well characterized clean Cu(OOl) surface; which will also serve as a 
reference to appraise the extent to which CO chemisorption affects it. Notice in Table [T] 
that both LDA and GGA-PBE produce an inwards relaxation of the surface layer of ~ 3%. 
These results are in agreement with earlier pseudopotential calculations^i^>2Ii^i^ and with 
surface structure measurements via the medium ion energy scattering (MEIS) technique (see 
Table [H).^ 

As illustrated in Fig.[Tt^a), the primitive super-cell of the CO-covered surface contains two 
Cu atoms per layer, which are non-equivalent in odd numbered layers since CO sits directly 
above only one of them. Accordingly, atoms in the first layer are referred either as covered 
or bare atoms. To our knowledge, there is no experimental characterization of the substrate 
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geometry after CO adsorption. Nevertheless, in agreement with previous calculations,-*^ 
our results^i^S indicate that the first and third layers rumple, while the second layer atoms 
do not. Quantitative differences, however, arise between LDA and the GGA-PBE results 
regarding the interlayer relaxations of the Cu(OOl) surface upon CO adsorption. We refer the 



reader to Table 1 in Ref. 
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for the details of the structure of c(2x2)-CO/Cu(001) and their 
comparison with available experimental data and other calculations.->21 Here, we only stress 
that, although both LDA and CCA agree on the fact that CO raises the original inward 
contraction of covered atoms and even makes them relax slightly outwards with respect to 
the bulk situation, such outward relaxation is two times larger within the GGA-PBE than 
what it is within the LDA. In addition, according to LDA, the inward relaxation of bare 
atoms is slightly decreased by CO, whereas GGA-PBE predicts that bare atoms undergo an 
inward relaxation which is even larger (~ -4.0%) than that of the topmost atoms of the clean 
surface. Note in passing that for the Cu-C and C-0 bond lengths, for which experimental 
data is available,— GGA-PBE gives slightly better agreement to the experimental results 
than LDA does.^^'^° Likewise, our GGA-PBE calculation produces a chemisorption energy 
of 0.68 eV at ^ = 0.5 ML, to be compared to the experimental valued of 0.57 eV. LDA, on 
the other hand, gives a much higher value of 4.00 eV. 

IV. RESULTS AND DISCUSSION OF THE LATTICE DYNAMICS 

Just for reference, we mention that our calculated phonon dispersion of bulk Cu 



(see Ref. |20| ) is in reasonable agreement with neutron inelastic-scattering measurements 
(NIS).— 1^ As for the dispersion of the surface phonons of Cu(OOl), they have been stud- 
ied at length recently by DFPT methods using both LDA and GGA-PBE.— Neverthe- 
less, we have repeated such calculations for consistency in comparison with results for the 
chemisorbed surface, which is the main subject of this work. The notation used in this paper 
for the surface phonons is in accordance with that introduced by Allen et al.^- The (1x1) 
and c(2x2) SBZ corresponding to the clean and the chemisorbed surfaces, respectively, are 
shown in Fig. [1] (b). We reiterate that the (1x1) SBZ is defined so that its zone boundaries 
along the [100] and [110] directions correspond to the M and X points, respectively. In 
addition, the T-M, T-X, and X-M segments are denoted as the S, A and Y directions, 
respectively. Before proceeding with our analysis, it is important to notice that the q-space 



8 



of (1x1) SBZ which is outside the c(2x2) SBZ is back-folded within the c(2x2) SBZ. 

We now turn briefly to the highhghts of the results of our calculations of the phonons of 
Cu(OOl). Full details of our results can be found in Ref]20l. Note that although a number of 
surface modes and resonances are found — especially along Y —,22, we describe below only 
those modes which are of interest in the discussion of the chemisorbed surface (see Figs.[2]and 
[3]). Table HTl summarizes the Cu(OOl) phonon frequencies at high symmetry q-points {X and 
M) and compares with those reported in Ref. Il3| and experiment.— We should point out 



that the sharp longitudinal resonance detected in HAS measurements,— is not reproduced 
in DFPT calculations. An explanation as to why theory and experiment differ on this 
particular issue remains an open question that deserves further investigation. Nevertheless, 
for the purpose of this work, we shall not pursue the origin and/or modifications of such 
longitudinal resonances since they completely vanishe after CO chemisorption.— Notice that 
the issue of the longitudinal resonance persists for several metal surfaces.— Chis et al. have 
recently addressed the case of Al(OOl)^ and Cu(lll).— 

With regard to the c(2x2)-CO/Cu(001) system, as discussed above, LDA indicates that 
the chemisorption of CO has little impact on the force constants mediating the interaction 
between hare atoms and their first NN in the second layer, whereas those mediating the 
interaction between covered atoms and their first NN in the second layer, YY, YZ, ZY, 
and ZZ, are mildly softened by 12, 12, 24 and 6%, respectively, as compared to the clean 
surface.— GGA-PBE, in contrast, finds that CO chemisorption modifies the force constants 
of not only covered atoms but also of hare ones. Naturally, the major effect occurs on 
covered atoms, whose corresponding force constants are softened by 40, 20, 38, and 14%, 
respectively, while those of hare atoms are stiffened by 14, 14, 14, and 9%, respectively. 
Unlike the interlayer force constants, intralayer force constants between hare atoms and 
their first NN covered atoms are barely altered by CO chemisorption. Namely, they are 
~5%(GGA-PBE) or ~8%(LDA) stiffer than would be if CO were not present. Tables [m] 
and IIVI summarize the frequencies obtained at F and X, respectively, and compare them 
with those found in experiments, when available, and former theoretical studies. In the 
sections below, we discuss in detail the characterization of the vibrational modes displayed 
by the Cu(OOl) substrate and shown in Fig. [31 Note, however, that because of the failures 
of LDA outlined in Section I and for the sake of clarity, our discussion will be focused on 
the results obtained via GGA-PBE. LDA results are mentioned only for cases of substantial 
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discrepancy with GGA-PBE. 
1. Substrate modes along S 

In the following analysis, one must bear in mind that modes proper to the clean Cu(OOl) 
surface (i.e. within the (1x1) SBZ) along the S direction (from the F point to the M point) 
which appear in the latter half of the same are now accessible in the first half, which lies 
inside the c(2x2) SBZ (see Fig. [I]). 

(a) The Si mode. The Rayleigh wave, known as 5*1,^^ is the surface mode with the lowest 
frequency along S direction. On the clean surface, at M, its polarization is mainly vertical 
(V) and localized in the first-layer, acquiring an additional longitudinal (L) polarization as 
it approaches the F-point. On the chemisorbed surface, 5*1 increases its frequency along S, 
until it reaches the zone boundary and matches, except for a 1 meV gap, the back-folded 
Si, which has its maximum at F and decreases its frequency along S. Regarding the gap 
between the RW branches at the zone boundary, we see that the higher (lower) branch 
corresponds to a mode whose amplitude weight is primarily vertical in bare atoms {covered 
atoms dragging the molecule to some extent), involving also a small contribution from the 
L-polarized vibration of covered (bare) atoms. In fact, the higher branch matches the RW of 
the clean surface at the zone boundary, suggesting that bare atoms are not affected by the 
presence of CO. Nevertheless, we will see that this is not always the case. At F, where the 
softening with respect to the clean surface is maximal, the back-folded Si mode corresponds 
to an out-of-phase vibration between covered and bare atoms. In this case, the contribution 
of bare atoms is 50% larger than that of covered atoms and CO molecules are dragged 
parallel to the latter. Close to F, Si is broadened and appears as a finite-width resonance 
whose maximum amplitude weight can be as low as 6% in the first two layers. It softens by 
~16% at F, overestimating in fact HAS measurements (~10%). We observe that although 
LDA and GGA-PBE predict different effects of CO chemisorption on the force constants of 
the first layer, in both cases 5*1 softens by the same percentage. 

(b) The S2 mode. This mode is mostly V-polarized and localized in the second layer 
at the zone boundary. On the clean surface, this mode soon forms part of a band of bulk 
resonances along S whose maximum amplitude weight comes from either the second and 
first layers (with V- and L- polarization, respectively) when close to the SBZ boundary, or 
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first layer (with mixed V-L polarization) when close to F. GGA-PBE and LDA yield slightly 
different results as to its prevalence along S and its degree of localization (see RefQ. 

On the chemisorbed surface, however, 5*2 only appears back-folded from M to F. It steeply 
disappears as soon as it immerses into the bulk band. At F, 5*2 is well inside the bulk band, 
even though, no coupling to bulk modes is observed. On the contrary, it is more localized 
on the chemisorbed surface than on the clean surface (at M). S2, in addition, stiffens by 
~ 3 meV on the chemisorbed surface. 

(c) The Li mode and the corresponding shear-horizontal (SH) branch. Li refers to a L- 
polarized mode mostly localized in the first layer of the clean surface at the zone boundary 
of the SBZ. It rapidly becomes a resonance as it approaches F. As a consequence, on the 
CO-covered surface, it appears mainly as a back-folded mode. At F, the back-folded Li 
and its degenerate SH pair stiffen very slightly (~ 1%) due to the chemisorption. Neither 
Li nor its degenerate SH pair are totally localized in the first or second layers (as in the 
clean surface) but exhibit a significant contribution of ~60% to the amplitude weight from 
deeper layers. Outside F, the back-folded Li rapidly disappears, while the back-folded SH 
mode disperses and becomes more localized towards the zone boundary, where it matches 
the branch that originally appears on the clean surface but vanishes rapidly back to F. 



2. Substrate modes along A 

Here, one should notice that, like the back-folding along S, modes proper to the clean 
surface along Y are back-folded to A in the c(2x2) SBZ (see Fig. [U (b)). 

(a) Si and the corresponding L-branch. Si is the mode of lowest energy of Cu(OOl) and is 
totally localized in the first layer at the X point. It cannot be detected by standard planar 
scattering techniques on the clean surface since its polarization is SH all along A. 5*1 changes 
rapidly to V- polarization along Y but, to our knowledge, no experimental data along Y is 
yet available. 

On the chemisorbed surface. Si is totally localized in the first layer and the CO overlayer. 
The latter vibrates in-phase with the covered Cu atoms but with much smaller amplitude. 
Si softens at X by ~8% (see Tables HI] and HV]) . Interestingly, the calculations reveal a 
L-polarized mode that is degenerate with Si at X, albeit the degeneracy is broken outside 
X. Such mode corresponds to the section of Si along Y in the SBZ of the clean surface. 
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which is back-folded along A on the chemisorbed SBZ and changes polarization from SH to 
L. The back-folded Si remains L-polarized as it goes across A towards T up to the crossing 
point with the RW, where it becomes V-polarized. It is slightly softened around the zone 
boundary, just as much as Si at X; nonetheless, the softening becomes stronger — similar 
to that of the RW — right after crossing the RW and the transitioning to V-polarization. 
Note that Ellis et al. observed some HAS peaks precisely at the region where back-folded 
Si is V-polarized.— Those peaks were at the time said to be associated with multi-phonon 
excitation bands. The excellent fit of their measured dispersion to our back-folded Si {Y 
is back-folded onto A), however, suggests that this latter mode is observed rather than 
multi-phonon excitation bands. 

(b) The 5*4 mode. This is the RW along A. On the clean surface, it is essentially V- 
polarized, although a L-contribution is also present. Although its amplitude weight is great- 
est at the first layer, it decays slowly as a function of the layer depth. On the chemisorbed 
surface, it is also mostly localized (~60-70%) in the first layer and, to lesser degree, in 
the molecule. 5*4 also displays a small splitting at X (see Table HV]) . yet in this case both 
branches are softened with respect to that of the clean surface by 9.6 and 15.0% at X, 
while HAS measurements^^ find a softening of ~8.2%. From the first layer, only the covered 
(bare) atoms contribute to the mode with lower (higher) energy. The back-folded S4, on the 
other hand, originally arises along Y on the clean surface and has V-polarization at the zone 
boundary and a predominantly SH-character as it crosses Y. On the chemisorbed surface, 
this mode is also V-polarized close to the zone boundary. In GGA-PBE, it broadens and 
becomes a resonance as soon as it immerses into the bulk band, reappearing as a surface 
mode close to F with mixed L- and SH-polarization. In LDA, however, S4 remains highly 
localized on the surface and changes smoothly to L-polarization.— 

(c) The 5*5 mode. The polarization of this mode is SH and it is predominantly localized in 
the second layer on the clean surface. On the chemisorbed surface, it is also strongly localized 
in the second layer but with mixed SH- and L-polarization. It rapidly becomes a resonance 
along A towards F. At X, 5*5 splits (see Table \W\ . In GGA-PBE both branches are totally 
localized in the second layer. One of these softens by 8.4% and the other stiffens by 3.5%. 
LDA shows that the lower branch softens while the other - bearing ~25% contribution from 
deeper layers - does not change at all. 

(d) The Sq mode. This mode developes inside the spectrum gap, close to the zone 
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boundary with a frequency of 23.6 meV (see Table II V|) . On the clean surface, it has a 
predominant L-polarization in the first layer. On the chemisorbed surface, however, it 
is found to be degenerate with a SH pair at X. Note in Fig. [3] that the degeneracy is 
broken inside the SBZ. Our calculations find that 5*6 - and the SH-branch - remain nearly 
dispersionless and involve an additional V polarized second layer vibration, as found on 
the clean surface. These modes are, incidentally, more localized in the first layer of the 
chemisorbed surface than in that of clean Cu(OOl). At the right end of the bulk band gap 
(see Fig. [3]), 5*6 slightly softens and becomes a resonance while it penetrates the bulk band. 
The SH branch, in contrast, extends well inside the bulk band and becomes more localized 
at the top two layers. According to GGA-PBE, 5*6 softens at X by 2.0%, while LDA predicts 
no softening (see Table HV]). 



V. SUMMARY 



A first-principles study of the dynamics of a c(2x2)CO overlayer chemisorbed on Cu(lOO) 
has been presented. Our calculations show that LDA displays a mild effect of CO chemisorp- 
tion on the relaxations and force constants of the surface. The LDA predicted adsorption 
site is also not in agreement with experimental results. Nevertherless, LDA is able to give 
good agreement with the HAS data for the frequency at the zone center of the hack-folded 
RW. This feature is rather unexpected in consideration of the poor ability of LDA to de- 
scribe the Cu-CO interaction and its less successful description of the acoustic modes of bulk 
Cu and the RW in Cu(OOl).— GGA-PBE, on the other hand, finds that CO chemisorption 
significantly perturbs the structure and the first NN force constants of the surface layer 
atoms. Softening of the RW is well reproduced; only slightly overestimated by ~1 meV at 
F. It is surprising that, while LDA and GGA-PBE differ considerably in the response of the 
substrate to CO chemisorption, the actual percentage softening of the hack-folded RW with 
respect to the clean surface is comparable. We find that softening of the RW along S and 
A is not necessarily connected with the vibration of the covered atoms, which indicates that 
mass overloading alone cannot account for it. If that were to be the case, all surface modes 
involving the first layer atoms would be softer. On the contrary, Li, for example, does not 
soften, notwithstanding the leading involvement of first layer atoms. In fact, it slightly stiff- 
ens, consistent with the hardening of intralayer force constants of the first layer. Moreover, 
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softening/stiffening of the various modes cannot be simply dictated by the atomic layer, 
propagation direction, and polarization characteristics, indicative of complexity in the mod- 
ifications in Cu force constants. For example, 5*2, the V-mode in the second layer, stiffens. 
In turn, while the 5*1 mode (along A) slightly softens in the region where it is L-polarized 
and totally localized in the first layer, it undergoes a stronger softening in the region where 
it is V-polarized and involves contributions from deeper layers. Softening of the RW seems 
thus due not only to the mass of CO but also to interactions of longer range, i.e. beyond 
first NN and involving deeper layers, that subdue differences in the local bonding of surface 
atoms and result in an overall softening of the RW, which is sometimes independent on 
whether covered or bare atoms are primarily involved. 

Our results call attention to the importance of the folding of Y (of the 1x1 SBZ) onto A 
(of the c(2x2) SBZ), which displays back-folded modes along the latter. For example, back- 
folded Si and back-folded S4 are found along A with changed polarization that may make 
them observable by standard scattering techniques. In particular, 5*1 acquires V-polarization 
close to r and L-polarization close to X. The good agreement for the dispersion of this latter 
mode HAS measurements^^ is convincing evidence that these peaks arise from the back-folded 
Si mode, rather than from multi-phonon excitation bands. Perhaps even more importantly, 
back-folded Si becomes as well discernible to planar scattering spectroscopy techniques close 
to the zone boundary since its polarization changes from SH- (along Y of the (1x1) SBZ) to 
L- polarization (along A of the c(2x2) SBZ) in that region. The above thus implies that the 
frequency of 5*1 (SH-polarized in Cu(OOl)) can indirectly be measured at X via back-folded 
Si. We believe, in fact, that some of the peaks observed and assigned to 2T overtones^ 
peaks by Ellis et al. along A may instead correspond to back-folded Si, in the region where 
is L-polarized. Our results here thus call for a new interpretation of the HAS data. 
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TABLE I: Percentage change of the interlayer spacing of the outermost layers of Cu(lOO) compared 
to the bulk situation. 
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-2.4 
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45 " Ref. 33 



TABLE II: Frequencies (in meV) of the surface modes of Cu(lOO) at the high symmetry points X 
and M (see Fig. [1] (b)). The main polarization (SH, L, or V) is denoted in parenthesis and the 
superscript indicates the layer showing largest amplitude weight. 
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TABLE III: Frequencies (in meV) at T of the surface vibrational modes of c(2x2)-CO/Cu(001). 
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TABLE IV: Erequencies (in meV) at X of the surface vibrational modes of c(2x2)-CO/Cu(001). 
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FIG. [U (a) The top view of the surface shows CO (grey circles), and first (filled circles) 
and second (open circles) layer atoms of Cu(lOO). The 1x1 (dashed line) and the c(2x2) 
(solid line) surface unit cells are underhned. (b) The corresponding (1x1) (dotted line) and 
c(2x2) surface Brillouin zones (solid line) showing the F, X, and M points; and the A, S, 
and Y directions. 

FIG. |2) GGA-PBE phonon dispersion of Cu(lOO), modelled by a 50-layer slab. Theoret- 
ical surface modes (filled circles) are compared with HAS (open triangles) and EELS (open 



40 



41 



46 



circles) measurements taken from Refs. 

FIG.E]: GGA-PBE phonon dispersion of c(2x2)-CO/Cu(100), modelled by a Cu 50-layer 
slab. Note that the dispersion of the high-lying modes, which correspond to the C-0 stretch 
and the Cu-CO stretch modes, are omitted in this figure since the emphasis here is on 
the modes of the substrate. Filled circles denote theoretical surface modes. Experimental 



data are taken from Ref. Il6|: Filled circles and triangles were associated with multiphonon 
processes. Open circles correspond to the substrate Rayleigh wave. Squares were associated 
with the FT mode of CO on the perfect c(2x2) structure (filled) and on defects in the adlayer 
at lower coverage (open). 
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FIG. 3: 
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